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Absorption and emission in inorganic bipolar solar cells based on low dimensional structures 
exhibiting the effects of quantum confinement is investigated in the framework of a comprehensive 
microscopic theory of the optical and electronic degrees of freedom of the photovoltaic system. In 
a quantum-statistical treatment based on non-equilibrium Green's functions, the optical transition 
rates are related to the conservation of electronic currents, providing a quantum version of the 
balance equations describing the operation of a photovoltaic device. The generalized Planck law 
used for the determination of emission from an excited semiconductor in quasi-equilibrium is replaced 
by an expression of extended validity, where no assumptions on the distribution of electrons and 
photons are made. The theory is illustrated by the numerical simulation of single quantum well 
diodes at the radiative limit. 
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I. INTRODUCTION 



Knowledge about the luminescence of solar cells is of 
great value both for the characterization of the cell and 
for the analysis of the device performance due to its cor- 
respondence to the radiative losses. Established theories 
predicting the luminescent emission of bulk solar cells 
are mainly based on the generalized Planck and Kirch- 
hoff formulas [1^ relating the emitted photon flux to the 
absorptivity and the bias voltage. However, the gener- 
alized Planck spectrum was derived based on the cru- 
cial assumption that the carriers are well described by 
a local equilibrium distribution characterized by a sin- 
gle quasi-Fermi level and at the temperature of the lat- 
tice. This assumption is no longer justified in nanos- 
tructured solar cells based on quantum wells or quantum 
dots |2l3l4l5l6l7l8l9ll0lllj . The purpose of the present 
work thus is to provide a more general approach to ab- 
sorption and emission in bipolar semiconductor nanos- 
tructures where the critical assumptions are relaxed. To 
meet the requirements formulated above, a description of 
emission is derived within a general microscopic nonequi- 
librium theory of quantum-confined photovoltaic systems 
[12]. The appeal of the approach lies in the way the op- 
tical properties of a semiconductor nanostructure are re- 
lated to its electronic properties on very general grounds. 
Based on non-equilibrium quantum-statistical mechanics 
and starting from a quantum theory of electrons and 
photons in semiconductors nanostructures, absorption 
and emission are calculated from the transverse polar- 
ization function, which appears as the convolution of 
the non-equilibrium Green's functions for electrons and 
holes [1 3I14I15I16I17I18I19I2QI2T] . Unlike more basic ap- 
proaches, the theory allows the inclusion of many-particle 
corrections and interaction effects, such as those due to 
electron-phonon coupling or exciton formation. 



II. THEORY 

The following outline of the general theory of quantum 
photovoltaic devices is restricted to the radiative limit, 
i.e. non-radiative recombination (Shockley-Read-Hall, 
Auger) and generation (impact ionization) processes are 
not considered. 



A. Hamiltonian 

At the radiative limit, the full quantum photovoltaic 
device is described in terms of the model Hamiltonian 

n =He + + Hp, total (1) 

He =Hl + electronic (2) 

TYg =7Ye7 + + ^ee, interaction (3) 

consisting of the coupled systems of electrons (TYe), pho- 
tons (7^7 ) and phonons (TYp). Since the focus is on the 
electronic device characteristics, only He is considered 
here, however including all of the terms corresponding to 
coupling to the bosonic systems. Within the electronic 
part, TYg contains the kinetic energy, the (bulk) band 
structure and band offsets, and it also includes the elec- 
trostatic potential from the solution of Poisson's equa- 
tion. The interaction part H\ consists of the terms 7Ye7, 
Hep and Hee for the interactions of electrons with pho- 
tons, phonons and electrons or holes, respectively. While 
the first term describes generation and recombination of 
free charge carriers via the absorption and emission of 
photons, the second one provides the carrier relaxation 
via the absorption and emission of phonons. The last 
term contains the carrier-carrier interactions responsible 
for screening effects or exciton formation. 

Within the non-equilibrium Green's function theory 
of quantum optics and transport in excited semiconduc- 
tor nanostructures, physical quantities are expressed in 
terms of quantum statistical ensemble averages of single 
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particle operators, namely the fermion field operator ^ 
for the charge carriers, the quantized photon field vector 
potential A for the photons and the ionic displacement 
field U for the phonons. The corresponding Green's func- 
tions are 



G(12) 
^7.(12) 



:(*a)*^(2))t 



-{Ml)Au{2))c. 



electrons (4) 
photons (5) 
phonons (6) 



where (...)c denotes the contour ordered operator aver- 
age peculiar to non-equilibrium quantum statistical me- 
chanics [13 14J for arguments 1 = (ri,ti) with temporal 
components on the Keldysh contour [14]. 



B. Dyson's equations 

The Green's functions follow as the solutions to corre- 
sponding Dyson's equations |16|17|22j . 
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Go, Dq and Dq are the propagators for noninteracting 
electrons, photons and phonons, respectively, ^ denotes 
tensorial transverse quantities and j labels the phonon 
mode. The electronic self-energy S encodes the renor- 
malization of the charge carrier Green's functions due to 
the interactions with photons and phonons, i.e. gener- 
ation, recombination and relaxation processes. Charge 
injection and absorption at contacts is considered via an 
additional boundary self-energy term reflecting the open- 
ness of the system. The photon and phonon self-energies 
and describe the renormalization of the optical 
and vibrational modes, leading to phenomena such as 
photon recycling or the phonon bottleneck responsible 
for hot carrier effects. 



C. Microscopic optoelectronic conservation laws 

The macroscopic equation of motion for a photovoltaic 
system is the continuity equation for charge carriers, 

StPc(r) + V•jc(r) = ac(r)-7^c(r) c = e,h (10) 

where pc and jc are particle and current density respec- 
tively, Qc the generation rate and IZc the recombination 
rate of carriers species c. The corresponding microscopic 



conservation law reads [13] 

Irni { {zh {dt, + dt,) - [Hoirr) - Ho{r2)]) G(12)} 

= lim / d3pe^(13)G(32) -G(13)Se^(32)] (H) 
Jc 



In steady state, the continuity equation (10) is reduced 
to 



V•jc(r) = ac(r)-7^c(r). 



(12) 



In the microscopic theory, the divergence of the current 
corresponds to the "<"-component of the RHS in (11) 



V-j(r) = -2e 



J 2nnJ 



Efjr,r';E)G<{r',r;E) 



S< (r, r'; E)G^{v', r; E) - G^(r, r'; E) 



X S< (r', r; E) - G<(r, r'; E)E^Jr' , r; E) 



(13) 



If the energy integration is restricted to one carrier 
species, the above equation corresponds to the micro- 
scopic version of (12), and expresses the total local ra- 



diative interband transition rate. The total radiative in- 
terband current is found by integrating the divergence 
over the absorbing/emitting volume, and is equivalent to 
the total global transition rate, and, via the Gauss theo- 
rem, to the difference of the current at the boundaries of 
the absorbing/emitting region. On the other hand, the 
global radiative rate equals the total optical rate R^^^, 
related to the global electromagnetic energy dissipation 
w through 



w 



2^ 



(14) 



and which, via the quantum statistical average of the 
Poynting vector operator, can be expressed in terms of 
photon Green's functions and self-energies [ 20f21) . pro- 
viding the following microscopic expressions for absorbed 
and emitted photon flux 

R'^P\t^)=R"^%{h^)-R';%^{t^), (15) 
K%sif^) = j '^^^ j [n>(r,r',fi^) -n<(r,r',fia>)] 
xD<(r,r',M, (16) 
i?°^*^(M= j d^'r J d^r'n<{r,r',hio)D<{r,r',hw), 



where the photon density of states 

= i{D> - D<) 



(17) 



(18) 



was introduced. The expression for R^^l^n replaces the 
generalized Kirchhoff formula [Ij for the photon flux 
emitted from an excited semiconductor. It allows for the 
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inclusion of a realistic electronic dispersion via the po- 
larization function 11^, reflecting the effects of electron- 
electron, electron-hole and electron-phonon interactions 
and the non-equilibrium occupation of these states, as 
well as the consideration of the optical modes specific to 
the geometry of the system. 

III. IMPLEMENTATION 

The general theory outlined above is implemented for 
quantum well layers inserted in the intrinsic region of a 
p-i-n diode. Details of this implementation can be found 
in [12]. As a first approximation, the charge carriers are 
treated as a system interacting with a specified equilib- 
rium environment of photons and phonons, i.e. only the 
electronic Green's functions are renormalized. The main 
aspects of the model for charge carriers, photons and 
phonons are given below. 

A. Charge carriers 

Electrons and holes are modelled within the framework 
of empirical tight-binding theory for layered semiconduc- 
tor heterostructures, where the carrier field operators are 
expanded in planar orbitals representing linear combina- 
tions of Bloch sums of localized atomic orbitals over the 
plane of periodicity. Interactions among carriers are con- 
sidered only on the Hartree level, equivalent to the so- 
lution of the macroscopic Poisson equation for the effec- 
tive potential. Inclusion of screening and excitonic effects 
would require the calculation of two-particle quantities in 
higher order perturbation theory. 

Carrier selective contacts are obtained by applying 
closed system boundary conditions at minority carrier 
contacts, which results in a complete supression of leak- 
age currents. 

B. Photons 

Light is described via free field modes with occupation 
determined by the combination of the spectra of the ex- 
ternal illumination and the black-body thermal equilib- 
rium with the environment. For practical computations, 
the absorption and stimulated emission due to the lat- 
ter can be neglected. The self-energy in the electronic 
Dyson equation describing the coupling to charge carri- 
ers is computed within perturbation theory on the level 
of a self-consistent first Born approximation (SCBA) to 
second order in the electron-photon interaction. 

C. Phonons 

Lattice vibrations are implemented within the har- 
monic approximation and for noninteracting acoustic as 



well as polar optical modes. The interaction with the 
electronic system is on the level of a coupling to an equi- 
librium heat bath and is described by the corresponding 
SCBA self-energy. 

D. Absorption, emission and photon flux balance 

The linear absorption coefficient a is related to the 
transverse polarization function or photon self-energy via 

a(M = -^^Ei"^n'pp(q = 0'M, (19) 

Im = -^( - U<), (20) 

where Sb is the background dielectric constant. For dipole 
coupling to free field photons, the self-energy reads 

nf(q;i?) / ^Tr{M^(k||,q,A)G§(k||;£') 

k|| 

xM^(k||,-q,A)G^(k||;E'-^)}, (21) 

where M^ describes the dipolar matter-field coupling and 
the trace is over both layer and orbital indices. For the 
considered model, the total radiative rate is given by 

Rrad = A ^ {</>I,q[n^(q. ^^q) - n^(q. ^^q)] 

q,A 

-n<(q,nc^q)}, (22) 

where A is the model layer width. From this representa- 
tion of the rate follows the intuitive result that the ab- 
sorbed flux is proportional to the incoming photon flux 
(j)^ , while the emitted flux only depends on the polariza- 
tion n<. 



IV. NUMERICAL RESULTS 

A. Models system: geometry and simulation 
parameters 

The simulated system consists of an ultrathin (50 nm) 
p-i-n diode with a single quantum well embedded in the 
i-region. A two band tight-binding model is used for 
the electronic structure with the the two bulk band gaps 
given by Eg^iow=0.6 eV and Eg^high=0.9 eV, with a sym- 
metric offset of A£^=0.15 eV. The effective masses and 
the parameters for the electron phonon-interaction are 
chosen for a GaAs/AlGaAs material system. The sys- 
tem is illuminated at a photon energy of 0.7 eV and an 
intensity of 1 kW/m^. 

B. Local density of states, absorption and emission 

Insertion of a quantum well in the z-region of a p-i-n 
diode locally extends the density of states below the high 
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FIG. 1: Local density of states for single quantum well p-i- 
n diode at vanishing transverse momentum and zero applied 
bias voltage. The effective band gap in the absorbing region 
is determined by the separation of electron and hole confine- 
ment levels. The interference pattern of the LDOS reflects the 
closed system boundary conditions at minority carrier con- 
tacts. 
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FIG. 3: The absorption exhibits a significant contribution 
from the quantum well. There is a pronounced blue shift in 
the fundamental emission edge with increasing voltage due to 
decreasing field (Franz-Keldysh effect), with corresponding 
increase of the effective band gap as observed in the local 
density of states. 



band gap bulk edge, red-shifting the effective absorption 
edge, which is now given by the position of the lowest con- 
finement level, as can be seen in Fig. [T] showing the local 
density of states at zero bias. Under very high fields, as 
in the case considered, the quantum well states hybridize 
with the bulk continuum outside the well. The result- 
ing significant deviation from the 2D-D0S is reflected in 
the absorption, shown in Fig. [3j This so called Franz- 
Keldysh effect includes an additional red-shift in the ef- 
fective band edge, which decreases with the effective field 
at the quantum well location under the application of an 
external bias voltage and resulting in a blue-shift of the 
emission edge as the voltage increases from V = V to 
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FIG. 2: Local density of states a large bias voltage close to 
Voc- The effective band gap exceeds the zero bias value due 
to a less pronounced Stark effect at the lower effective field. 



V = 0.44 V ^ Voc, which is displayed in Fig. [2] 



C. Local radiative interband rates 

Fig. |4] shows the local interband rates for genera- 
tion, recombination and the sum of the two effects, as 
computed via the electronic Green's functions and self- 
energies, consisting of inscattering (=generation) and 
outscattering (=recombination) components. At illumi- 
nation photon energies exceeding the effective band gap 
but below the higher bulk gap, net generation and re- 
combination are spectrally well separated and spatially 
confined to the quantum well, i.e. the region of finite 
overlap of absorbing electronic and hole quantum well 
states. On the other hand, integration over energy pro- 
vides net interband current of the entire device. 



D. Local current spectrum 

The spectral difference in generation and recombina- 
tion rates leads to a spatial inhomogeneity in the gen- 
eration and recombination currents, as seen in Fig. [5j 
with localization of the maximum in the quantum well, 
even though the integrated current is locally conserved, 
and the separate contributions of radiative dark and pho- 
tocurrent grow towards the contacts and reach their max- 
imum outside the absorbing region. 
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preaches for the modehing of absorption and lumines- 
Quantum effects in low dimensional nanostructured cence. A comprehensive microscopic theory of quan- 
absorbers for high efficiency solar cells require new ap- tum photovoltaic devices including the effects of quan- 
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turn confinement on the optical, electronic and vibra- 
tional properties can be formulated based on the non- 
equilibrium Green's functions theory of excited semicon- 
ductor nanostructures. Within this theory, microscopic 
electronic and optical conservation laws are derived and 
interrelated. Numerical simulations based on an imple- 



mentation of the theory for quantum well pin-diodes pro- 
vide a picture resolved in space and energy of the gen- 
eration, recombination and the corresponding currents 
in quantum well absorbers at the radiative limit and for 
both short circuit and close to open circuit conditions. 
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